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Abstract 

The control of systems witli unknown dynamics and unpredictable disturbances has raised some challenging 
problems. 'This is particularly important when high system peformance is to be guaranteed at ail times. Recently, 
the Time Delay Control has been suggested as an alternative control scheme. The proposed control system does 
not require an explicit plant model nor does it depend on the estimation of specific plant parameters. Rather, it 
combines adaptation with past observations to directly estimate the effect of the plant dynamics. 

This paper formulates a control law for a class of dynamic systems and then presents a sufficient conditiorf for 
control system stability. The derivation is based on the bounded input-branded output stability approach using 
Loo function norms. The control scheme is implemented on a five degrees * f- freedom high speed and high precision 
magnetic bearing. The control performance is evaluated using step respond s, frequency responses and disturbance 
rejection properties. The experimental data show an excellent control perf >i inance despite t he system complexity. 

1 Introduction 

Some classical control methods deal with well known linear time-invariant systems. In many applications, however, 
some relevant part of the system maybe unknown, time varying, or nonlinear. Controlled systems are thus often 
limited to operating in only a small portion of their available range. For example, servo motors must operate in 
the linear part of their range for accurate control. Restrictions such as these have led to the developrn* nt of control 
techniques that deal with such complexities. 

Several types of modern control strategies have been developed to deal with nonlinear, time- variant systems. One 
of the first methods to accommodate nonlinear systems was Model Reference Control. This technique employs a 
model of the system and uses the difference between the model response and the plant response as the input signal to 
the plant [18]. The model is either a physical model or a simulated system on a computer. Although it has no variable 
parameters, it is very useful for either specifying desired performance or for the observation of unaccessabie states. 
A drawback in this technique is that it requires knowledge of the full dynamic model and system limitations. When 
perfect cancellation of the system nonlinearities is not achieved due to imperfect modeling or inaccurate parameter 
values, the dynamic performance of the plant may be degraded to the point of closed loop instability [22]. 

Another advanced technique is Adaptive Control. An adaptive system measures a certain index of performance 
which is a function of the inputs, states and/or outputs of the system. From the comparison of the measured index of 
performance with a set of given ones, the adaptation mechanism modifies the parameters of the controller or the set of 
given ones [4,14,16,21]. 1 here are several classes of adaptive control. A very common variation uses a desired reference 
model as a basis for comparison and is termed Model Reference Adaptive Control (MRAC). In the direct MRAC, no 
attempt is made to identify the plant parameters. Controller parameters are directly updated. In Self- Tuning control, 
plant model parameters are identified/modified and the controller action is automatically updated according <>o a 
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2 Time Delay Control 

In this paper we are concerned with a class of systems described by the following differential equations, 


x(t) = F(x, 0 + G(x, u, t) + D(<) (1) 

where x(*) E 7v n and u(f) E are the system state vector and control input vector respectively. F(xJ), G(x,uJ) 
and D(0 are vector functions with appropriate dimensions and represent respectively known dynamics, unknown 
dynamics and disturbances. The variable t represents time. In order to transform the system into a familiar form, 
Equ.(l) can be written as 

x(0 = F(x, t) + H(x,u,*) + Bu(0 (2) 

where the new term H(x,u,f) is defined as 

H(x, u, t) = G(x, u, t) 4- D(*) - Buff) (3) 


and B is a matrix to be selected by the designer. A reference model that generates the desired trajectory is chosen 
as a linear time invariant system, 

x m (<) = A m x m (0 + B m r(0 (4) 

where x m (/.) E 7v n is a reference model state vector, r(t) E U r is a reference input. A m and B m are constant 
matrices with appropriate dimensions. 

The class of systems considered in this paper includes systems that satisfy a matching condition. It w'as shown 
m references [30,32,33] that systems in a special canonical form satisfy the matching condition. These systems can 
be partitioned as follows 
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where the partial states are x, 6 TC n_r ,x r € K r ,x a = [x r+ i,x r+2 . . . ,x„]' 6 K n r ■ The vector functions have the 
following dimensions F r (x,<), H r (x,u./) G H r ,B G tt nXr and B r € ft rxr is of rank r The matrices involved in 
the reference model of Equ. (4) are also partitioned in the same manner, 
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where K G lZ nxn and K r G K r * n 
according to 


The objective is to generate a control action u that forces the error to vanish 


e = (A m + K)e = A e e 


( 5 ) 


The control action that combines past observations with adaptation for systems described by Equ (2) is given by 


ii(/) = B + [— x(< - L) + F(x,t - L) - F(x, <) 
+A, n x(<) + B m r(4) + Bu(< - L) - Ko.(t)} 
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where the parameter L represents the time delay [30,32,33], the error vector e is defined as the difference between 
the plant and the reference model state vectors, 
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3 Stability Analysis 

3.1 Error dynamics 

A -a- . -i in m nl the stability of such control systems using time delay depends on the delay parameter L, the 
As indicated in the staouiiy oi / c _„j n c .p^nonse of the reference trajectory. The 

tt;rde‘5 in^'Sltpn. stability procedure .M 

follows, w. discus, the stability analysis for 
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+B r ii(< — L) + A mr x(/) + B mr r(<) - K r e(<)] 


H r (x,u,t) - H r (x,u,t - L) + A mr x(<) + B mr r(<) — K r e(0 J 

The previously defined error e of Equ.(7) is now governed by 

c (0 - (Am + K)e(!) + H(x,u,< — L) - H(x,u,<) (9) 

where the second and third tern, arc forcing functions due to the unknown system dynamic, and unpredictable 


disturbances. 


Rewriting Equ. (9) as 
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and 


e(|) = (A m + K)e(0 + p(0 


p (0 = 
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P r(t) 


( 10 ) 


p r (*) = H r (x,ud — L) — H r (x,u,0 

One may ask the question; what conditions does the vector p (t) have to satisfy for the system to be stable? A 
sufficient condition for stability will be derived in the next section. 

3.2 Sufficient conditions for stability 

This section presents a general solution to this multi-input multi-outp it control problem We will use the bounded 
input-hounded output approach based on L 0 0 norms in order to derr e sufficient conditions for stability. We now 
consider the governing differential of the error as given by Equ.(10) an I its correponding time response, 

e(<) = e< A ”* +K >‘e(0) + [‘ e (A - +Kl ‘- r| p(r)</7- 
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We will use ||( )r|| indicate the norm of the time trancated function (.) and ||(.)||, f° r the induced i latrix norm. 
Taking the norm of the error [28], 
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The desired error dynamics given by ( A m + K) are always chosen to be asymptotically stable. This implies that 
there exist finite positive constants m, A such that 
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Therefore, the norm of the error is bounded 


||e T || < a + /?||pr|| l 11 ) 

where 

a = m|| o(0)||, 0 = y 

In order to be more specific on these stability conditions, we need to expand the forcing term p(0 We can rewrite 
Pr(<) as 



Pr(0 = 


H r (x(< - L),u(t - L),t- L)-H r (x(t),u(t - L),t) 

+H r (x(t),u(t - L),t) - H r (x(t),u(<)./) 

||[p(<)]t|| = ||br(0]Tll<ll[Hr(x(t-I),u(«- L),t - L) (12) 

-H r (x(i).u(< - L),0 ]t|| 

+||[H r (x(t),u(< - L),t) - H r (x(t).«(0.0 rll 


Assuming that the function H r (x,u,<) is continuous and differentiable, the Mean Value Theorem yields [6] 
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Equations (13) and (14) involve terms in x and u 
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To express them in terms of e, the following expressions are 
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since e r = x mr — x r and x mr (0 — 
algebraic manipulation yields, 
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where i G (< - L,t). Substituting Equ. (18) in Equ. (17) yields 
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||Pt|| £ c 1 +C2||e7>|| +c 3 ||eT|| 
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where the constants ci,C 2 and c 3 can be identified as 
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We will assume later that the terms in the right hand side of the three equations listed above are bounded and hence 
the parameters <?i, c 2 and c 3 will be bounded. The norm of the error in Equ.(ll) can now be evaluated and is found 
to be, 


||er|| < « + 0 IIptII < a + f)ci+ /?o|| er || + 0c 3 || e T || 
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The above equation relates the norms of the error and the error derivative To obtain absolute bounds on the error 
another equation of this form is needed. The norm of the derivative of t lio error e, can be found from Equs. (10) 
and (10), 
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This condition can be stated as, 
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and substituting Equ.(20) one obtains 
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the following sufficient condition follows, 

If the functions F(xJ) y G(x,uj),B(t) are continuous and differentiable, and 

If the Jacob, an matrices vector functions ^ and 

If the eigenvalues of the matrix (A m + K) are in the left-half plane, and 
If the following stability condition is satisfied 
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Then the time delay controller is stable, and 
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The constants c\ , cn , a and (3 are , 
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Figure 1: Region of Stability in the ||e 7 *|| — ||ex || space 


Condition (ii) implies that the rate of change of the functions F r ,G r and D r with respect t,x and u are bounded. 
Condition (iii) implies that the desired error dynamics are chosen to be stable. Condition (iv) relates the time delay 
the rates of change of F r and G r with respect to x and u, and the desired error dynamics specified by ( A m + K). 
The bounds on the norms of the error, e and its derivative e are shown graphically in Figure 1. Equ. (20) yields 
a straight line as the boundary between stable and unstable regions. Similarly, Equ. (22) yields another straight 
line. The common region bounded by the two straight lines is the region of stability as shown. The absolute bounds 
l|e|| max and ||e|| mar shown in Figure 1 correspond to those given by Equ. 1 25). 

3.3 Special cases and discussions 

Condition (iv) stated in ihe proposition 1 can be rewritten in a more convenient and usable form. Tins simplified 
and convenient condition, however, is more conservative. Starting with condition (24) and after some algebraic 
manipulations we obtain, 

„»G. n -, ... , { l-||i&II.M» + ?IIA,n-t-K||,l } 

11 wB ' (TTliWiiTT) T m) 

< + f [||A mr + K r ||, 1 

, +||A m + K ||i(l + Hex’ll* L )| I 

The condition (26) implies a bound on the variation of relative to the controller gain matrix B r . The size 
of this bound is dependent on the delay time L, the norms of Hop’ll* a nd ||^ 5 l ||«i an d the desired characteristics 
of error dynamics given by (A, n H- K). The smaller the delay time is and the smaller the bounds on ||#||< and 

are, then the larger the allowable size of bound is on the range of relative to B r . 

In the case of first order SISO systems the vector functions reduce to scalar functions G r = g, B r = b and F r = /. 
Some interesting results are stated below for this class of systems. For first order SISO systems, we have 

|| c (A m +K)(l-r)||. _ | € a m +k)(t-r) | _ e (a m + *)(*- r ) 
m — 1 

A = — (rtm + k) 
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||A mr +Kr||, = ||Am + K||i = I «m+*l= ~(*n + U 

the stability condition (26) becomes 
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being always positive. For the case where = b the stability condition (28) 
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intuitive sense. From (28), it can be observed that for sufficiently small l y >, 

the following limiting case, 
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The result indicates that stability is maintained for a variation of 66% of £ with respect to b. 

When the control distribution matrix # is a constant and known, the controller gam matrix B, may be chosen 
such that « = B This enab'es ex^t ^cdl^on » 

: i Tlii;r i - J. in,,,.,, o„ — - 

:””d 7t “? it the tili., condition becansc of exact canceUat.on. 

4 Application: Control of a high speed and high precision magnetic 
bearing system 


The magnetic extern node, I. . 
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. D r (0 . 


(31) 
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Figure 2: A five-axis magnetic bearing 


Rotor Mass = 2.2 Kg 

Air gap for thrust bearing = 400 /im 

Air gap for radial bearing — 250 fin l 

Maximum current to bearings = 10 Amp. 

Maximum rotational speed = 45,000 RFM 


Table 1: Relevant system parameters 


where x q G ft* and x r G ft* represent displacements and velocities of the rotor with respect to the bearing stators 
respectively. The current inputs to the electromagnets are represented by the vector u G ft* The control objective 
is to levitate the rotor and maintain stability. Also, the control system must reject disturbances under spinning and 
nonspinning conditions of the rotor. 

This plant is multi-input, multi-output with all five degrees of freedom unstable open loop. Disturbances and 
coupling include forces due to gravity, magnetic actions, unbalance and gyroscopic effects. All of these effects will 
show up in the vector function G r . Note that since the magnetic force is proportional to the current squared and 
inverseley proportional to the gap distance squared the function G r (x,u) is a nonlinear function dependent on the 
state x and the control action u. 

The variation of the component of G r in the Z direction in terms of the gap x and control current u are shown 
in Figures 3, 4, and 5. Figure 3 shows that for a gap of 0.15 mm, the rotor acceleration corresponds to 30 m/sec 2 
and 110 m/sec 2 for control currents of 1 and 2 amps respectively. Also, for the same gap opening and with currents 
levels of 1 and 2 amps, ff changes from about 0.22 xlO 6 ^ to about 0.88 x10 6 j£t, and §£ changes from about 32 
m/amp-sec 2 to about 112 m/amp-sec 2 . It is clear that this particular device experiences drastic dynamic changes. 
Therefore, such dynamic information would be necessary if a conventional control system is used: otherwise the 
system performance may be acceptable only for some specific operating conditions. 

The control algorithm was implemented on a DSP chip as shown in Figure 6. In this experimental setup, we have 
the option of controlling the system using either a linear analog controller which resides in the compensation block, 
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„ . . . : n DSP board The position signal for the Time Delay controller 

“ * d /'* y ^ TP2 » dM TM .~a lh» signal is .hen «. So an analog to digital convene. 

” Ih'ongh the test po.nU TP2 »d/o, 1 FJ a S Unstable built in low pass filter where 

££££&£: fit«d“^ ““.^gnl, - sen. on. through the D/A convene, wh.ch h„ a low 

^^LiStrit^C^Ien*!, language. The santplin, frequency used “ 

me pru^icu , . __ a am * 7 n. ic * r Thf» cutoff freauency of the filter lor the position 

computation time for ^ sensor is already bandlimited to 1 KHz. The position signals 

signa was ep board and the control actions representing currents are sent out through the D/A 

are then obtained » “* “ fife, with a cutoff frequency set to 1.5 KHz Although a sampling rate 

and^dampin, ratio of .707. ^^^'^^k/i^blT^^the^p^UoTiOf^ro^^nro^Tin^OOf^T^Ojirn 
tracks the reference mo e response ver error between the desired and actual position trajectories is 

sttv: r r:±C3:r.brt7S 
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drive amplifiers (Auxi y P g . , d to check the disturbance rejection properties, we 

^ il^ain aro^O rad/se. The static stiffness is about 100 MN/m 

reference posit.on and actud rotor poMtoo J* ^ ^ and ^ jt b spinning at 10900 RPM 20100 

PPM ^oVS and tum RPM When the rotor is not spinning, the static stiffness is about 200 MN/m and 
the minimum stiffness is about 500 KN/m It is clear ' tha “ he 

ZZ up hating . lower stiffness. These data demonstrate that such a 
control scheme possesses excellent robustness properties. 


5 Conclusion 

The time d.la, controller algorithm u.«a paa. observation, for adapt™ i» c»trcdh-* 

and unpredictable disturl ances. The time de ay con ro aw is . , bounded input-bounded output 

nonlinear input ac t,o. 0. J* - • hve-degree-of-freedom 

mfgScZrm'g ^e^nToTpeHormauce' evaluated using step response, and disturbance rejection properties. .» 
shown to be excellent despite the complex nonhneanties in the system. 
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Figure 7: Experimental time response data of the thrust bearing 
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